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We consider a simple model of a Josephson junction phase qubit coupled to a solid-state nanoelec- 
tromechanical resonator. This and many related qubit-resonator models are analogous to an atom in 
an electromagnetic cavity. When the systems are weakly coupled and nearly resonant, the dynamics 
is accurately described by the rotating-wave approximation (RWA) or the Jaynes-Cummings model 
of quantum optics. However, the desire to develop faster quantum-information-processing protocols 
necessitates approximate, yet analytic descriptions that are valid for more strongly coupled qubit- 
resonator systems. Here we present a simple theoretical technique, using a basis of dressed states, 
to perturbatively account for the leading-order corrections to the RWA. By comparison with exact 
numerical results, we demonstrate that the method is accurate for moderately strong coupling, and 
provides a useful theoretical tool for describing fast quantum information processing. The method 
applies to any quantum two-level system linearly coupled to a harmonic oscillator or single-mode 
boson field. 

PACS numbers: 03.67.Lx, 85.25. Cp, 85.85.+J 



I. INTRODUCTION 

Josephson junctions have been shown to be effec- 
tive qubit elements for solid-state quantum computing 
architecturesi2i2di£iii Several proposals for multi-qubit 
coupling introduce electromagnetioZi2i2ii2iiiii2ii&i2ii£ or 
mechanicaliSiii resonators, or other oscillatorsji2ii2i2i to 
mediate interactions between the qubits. Such resonator- 
based coupling schemes have additional functionality 
resulting from the ability to tune the qubits relative 
to the resonator frequency, as well as to each other. 
These qubit-resonator systems are analogous to one 
or more tunable few-level atoms in an electromagnetic 
cavity, and the dynamics is often accurately described 
by the rotating-wave approximation (RWA) or Jaynes- 
Cummings model of quantum opticsiSi 

For a qubit with energy level spacing Ae coupled 
with strength g to a resonator with angular frequency 
ujq and quality factor Q, the RWA is valid when both 
|cj — Ae/h\ -C wo/Q and g -C Ae. However, the reso- 
nant Rabi frequency, which is proportional to <?, is then 
much smaller than the qubit frequency Ae/h. There- 
fore, restricting g to be in the simpler weak coupling 
regime leads to quantum information processing that is 
slower than necessary, allowing fewer operations to be 
performed during the available quantum coherence life- 
time. 

The threshold theorem 22 i 23 i 24 i 25 states that if the com- 
ponent failure probability p is below some threshold p t h, 
a computation with an error probability bounded by rj 
may be accomplished, provided a sufficient number of 
quantum gates are used for fault-tolerant encoding. In 
practice, it will be important to have p as small as pos- 
sible. To approach this limit, we wish to study qubit- 
resonator systems with stronger coupling (larger g) than 



may be correctly described by the RWA. This will allow 
us to consider faster switching times for qubit-resonator 
gates, and to understand to what extent the coupling 
may be increased while still retaining good fidelity. 

In this paper, we use a basis of dressed states^ to 
calculate the leading-order corrections to the RWA for a 
Josephson junction phase qubit coupled to a solid-state 
nanoelectromechanical resonator, or for any other model 
of a two-level system linearly coupled to a single-mode 
boson field. By comparison with exact numerical results, 
we demonstrate that the method is accurate for moder- 
ately strong coupling and provides a useful theoretical 
tool for describing fast quantum information processing. 

II. JUNCTION-RESONATOR DYNAMICS IN 
THE DRESSED-STATE BASIS 

A. Qubit-resonator Hamiltonian 

The Hamiltonian that describes the low-energy dy- 
namics of a single large-area, current-biased Josephson 
junction, coupled to a piezoelectric nanoelectromechani- 
cal disk resonator, can be written as^Siii 

H = 5Z £ ™ c m c m + ttooata - ig 2J x mm , c\ n c m , (a - a 1 *), 

m mm' 

(1) 

where the {cj„} and {c m } denote particle creation and 
annihilation operators for the Josephson junction states 
(m=0, 1, 2, . . . ), a and denote ladder operators for the 
phonon states of the resonator's dilatational (thickness 
oscillation) mode of frequency ujq, g is a coupling constant 
with dimensions of energy, and x mm > = (m\S\m ). The 
value of g depends on material properties and size of the 
resonator, and can be designed to achieve a wide range 
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FIG. 1: Two current-biased Josephson junctions (crossed 
boxes) coupled to a piezoelectric disc resonator. 



of values^ii An illustration showing two phase qubits 
coupled to the same resonator is given in Fig. ^ 

For simplicity we will consider only two levels in a sin- 
gle junction; generalization of our method to more than 
two junction states is cumbersome but straightforward^ 
However, all possible phonon-numbcr states are included. 
The Hamiltonian may then be written as the sum of two 
terms, H = Hjc + V . The first term, 



7tl 



c C\a' 



is the exactly solvable Jaynes-Cummings Hamiltonian, 
the eigenfunctions of which are known as dressed states. 
We will consider the second term, 

V = -ig[x 00 cic (a - a ] ) + x'oicjcia 

- x 01 c[c a^ + xuc\ci(a- a^)], (3) 

as a perturbation. The RWA applied to the Hamiltonian 
H amounts to neglecting V. Therefore, perturbatively 
including V is equivalent to perturbatively going beyond 
the RWA. 



where \mn) = |m)jj <8 |n)jj are the eigenstates of the 
uncoupled system. These states, together with |00), form 
a complete basis. The energies are 



W? = e + (j + l)^wo - — + o- — ^ — - ; 

and 77jc|00) = 0. On resonance, these reduce to 
\0,j + l}-ia\l,j) 



(8) 



V2 



(«d = 0) (9) 



and 



W° 



+ (j + l)huo + <ry/j + 1 



ffio(0) 



(wd = 0) 
(10) 

Below we will restrict ourselves exclusively to the reso- 
nant case. 

In what follows, we will need the matrix elements of V 
in the dressed-state basis, which are given by 



(2) mvw) 
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\J 3 + 2 ^oo S. 



3 + 1,3' 



3,3' + l 



+ a a' x n ( \/j + 16j+i,j' - V7 <5j,j'+i I (H) 



and 



igx 00 n (0) , 

{ipj \V\00) = — -=- 5 ja - cr dji. 



V2 



2V2 



(12) 



C. Dressed state propagator 



B. Dressed states 

It will be useful to define a set of Rabi frequencies 
according to 

fi> d ) = v /[%(0)] 2 +^ ; j = 0,l,2,... (4) 

where 

n J -(o) = 0' + i)*x2 fl |soilA (5) 

are the resonant Rabi frequencies for a qubit coupled to 
an oscillator containing 3 phonons, and where 

LUa = u! — Ae/h (6) 

is the resonator-qubit detuning frequency. The vacuum 
(i = 0) Rabi frequency on resonance is Oo(0) = 2g\xoi\/h. 

The eigenstates of -ffjc, or the dressed states, are la- 
beled by the nonnegative integers j = 0,1,2,... and a 
sign a = ±1, 

i^k = Pjj^jd) + ^d]|o, j + 1) - ianMlhj) (?) 

3 y^O,^) [fljiuJd) +o-uj d ] 



In quantum computing applications one will often be 
interested in calculating transition amplitudes of the 
form 

(f\e- im / h \i), (13) 

where |i) and |f) are arbitrary initial and final states of 
the uncoupled qubit-rcsonator system. Expanding |i) and 
|f) in the dressed-state basis reduces the time-evolution 
problem to that of calculating the quantity 

G${t) = (^\e- im ^:), (14) 

as well as (rf\e- lHt / h \00) and (00|e~ 1 00) . G$'(t) 
is a propagator in the dressed-state basis, and would be 
equal to 5 aa '5jjie^ iW i 1 *' ft if V were absent, that is, in 
the RWAiS Although it is possible to directly construct 
perturbative expressions for the propagator in the \mn) 
basis, the quantity defined in Eq. (|14f> turns out to be 
the simplest. 

To be specific, we imagine preparing the system at t = 
in the state |10), which corresponds to the qubit in the 
excited state m = 1 and the resonator in the ground state 
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n = 0. Wc then calculate the interaction-representation dressed states leads to 
probability amplitude 

c mn (t) ^ e^^( m n\e-^\l0) (15) <*>® = £ W0><00|e-^|^>, ^ 

for the system at a later time t to be in the state \mn). 

Here E mn = e m + nhuiQ. Inserting complete sets of the and. for mn =/= 00, 



Cmn (t) = e «W/»£; ^>->V f G % + f(^\!0) 



Using the relations 



and 



we obtain 



|0n> = ^[|VvU) + |C-i>] (forn^O) 



|ln) = -^[^+)-|V-)], 



c 01 (t) = l ^ Emtl 



it/>i AV ^ ^00 + G*dc 



00 

^00 ^00 



(17) 



(18) 



(19) 



(20) 



and 



cw{t) 



l e iE w t/h 



C++ c~\ — 

"00 "00 

— 1/ \G 00 G 00 



1 

-1 



So far everything is exact within the model defined in Eq. |JT]|, 

To proceed, we expand the dressed-state propagator in a basis of exact eigenstates \^ a ) of _ff , leading to 



(21) 



(22) 



Here £ a is the energy of stationary state |\I/ Q } . The propagator is an infinite sum of periodic functions of time. Wc 
approximate this quantity by evaluating the |^q) and £ a perturbatively in the dressed-state basis. 
The leading-order corrections to the dressed-state energies are of order V 2 . We obtain 



£ 



E 



L 00 



T 2 

x 01 



and 



m\vm\ 2 



w° 



v m\y\^)\ 2 

3=F3,° J J 



(23) 



(24) 



We will also need the second-order eigenfunctions, which, for a perturbation having no diagonal dressed-state matrix 
elements, are 



*oo) = A 



on 



^ (ib 1 ? \V\00) v (W\V\i>H,}{1>4,\V\00) 

i°°)-E^i^K)+ E 3 1 V r 

jct J jj'aa' 



(25) 
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and 



m + 



100) + 



(w:\v\^){ri"\v\ri 



E E 

j'< r ' 9 tjvj"< r "&<7 yVV 3 >^ V 3 



w° 



|00) 



(^/|y|oo)(oo|v|^; 
W?(W? - wp) 



(26) 



where A$n and the are normalization factors. 



Writing out Eq. (|2~2"|l explicitly as 

Gj/(t) = (VJ|*oo)(V^>oo)* e~ l£mtlh + E<^- l*5ff)(^<l*5ff)' e-^ t/fi , 



(27) 



and again making use of the fact that the matrix elements of V diagonal in j vanish, leads to 



5uu' A n „e 



-i£ 00 t/h 

+ ^00 ,™ W e 



+ A^(1-<W) 

+ 4l<(i-<w) 
+ E A 



(^>|oo)*(oo|v|^; 
«|y|oo)(oo|F|<) 



3° 



(^IV#><^o>l#>* 



(^'|i/|yf)*(^f|y|<)* 

_ (W '-Wr)(Wf-Wf) 

W|y|Vf)(VfM<) 
W-^W-wf) 



E 
E 



(Wf - Wg){W§ - W%') 



3^0,9 

-iej S t/n _|_ Q(y 3 ^ 



-i£otrt/K 



-i£ 0tr ,t/h 



(28) 



C aa '(t\ — ft , A 2 p-i£aat/h i a2 



2 2 

9 x oo 



2 W ff W ° 



-i£aat/h , ^2 

re -t- An 



1 



" 2 2 

9 ^oo 



x e 



-i£a„t/h 



^2 1 ^o-o-' 

0<T ' W CT - WK' 



o 



ff^oo y- / _J7 
2W n CT ' 4^ V 4(1 



0rT Wg - Wq' 

2 Y -Y ,_ 



2W n ° 



E 



9 9 

S &01 



g 2 X (7 g-X (7 i s 



(Wf - M^)(Wf - iy CT ') 



V4(W,f - Wf ) 2(W<f - Wf ) 



2 2 



2 2 

x oi 



VW-Wf) 2(W ff - Wf) 

-i£ 0tr 't/h 



(wi-wgyyvi-wf) 



-i£2o-t/h 



(29) 



where 



= V2 



(30) 



Note that there are no order V corrections to the dresscd-statc propagator. Because of this property, the leading order 
corrections are of order V 2 , and it is therefore necessary to use second-order perturbative eigenfunctions to obtain all 
such second-order terms. 

Finally, we note that the normalization constants are simply 



4>o = 



' oo 



1 01 



(31) 



A 0a = 



2 2 

9 x oo 



f(V2 



V (Wg - wfy 



2 x 01 



Wf) 2 



(32) 



n 2 r 2 
9 x oi 



2(Wf ) 2 ' 4 \ (Wf - W^') 2 



(V3 



^00 



V^crcr'sn) 5 



3 a; 01 



(Wf - Wf ) 2 



(33) 
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and 

A 2 a = 



2 x 01 



(VSxqo + V^crcr' xn) 2 (2 x 00 + a/3 aa' x n ) 



Ax 2 m 



(wf -wfy 



(ws - w^y 



(W 2 CT - Wf ) 5 



(34) 
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FIG. 2: Time evolution of probabilities |cio| 2 and |coi| 2 for 
the weakly coupled case of g/Ae — 0.03. Here the exact, 
RWA, and dressed-state perturbative results are essentially 
equivalent. 




FIG. 3: Probabilities |cio| 2 and |coi| 2 for the strongly cou- 
pled case g/Ae — 0.30. Here there are large deviations from 
the RWA behavior, which are correctly accounted for by the 
dressed-state perturbative method. Note the ten-fold increase 
in transfer speed compared with that of Fig. [5] 



III. TOWARDS INFORMATION PROCESSING 
WITH STRONG COUPLING 

In this section, we test our perturbed dressed-state 
method for the case of a finite-dimensional singlc-qubit, 
fivc-phonon system. The junction has parameters Ej = 
43.1 meV and E c = 53.4 neV, corresponding to that of 
Ref. 0. The resonator has a frequency ujq/2-k of 10 GHz, 
and the interaction strength g varies from weak (g <gC Ae) 
to strong (g f=a Ae) coupling. The bias current is chosen 
to make the the system exactly in resonance, and this bias 
is sufficiently smaller than the critical current so that the 
junction states are well approximated by harmonic oscil- 
lator eigenfunctions. The Hamiltonian for this system is 
diagonalized numerically, and the probability amplitudes 
Cmn (i) are calculated exactly, providing both a test of the 
accuracy of the analytic perturbative solutions and an es- 
timate of the range of interaction strengths g for which it 
is valid. Setting the initial state to be c mn (0) = S m i5 n o, 
as assumed previously, we simulate the transfer of a qubit 
from the Josephson junction to the resonator, by leaving 
the systems in resonance for half a vacuum Rabi period 

Figures |3 El and^J show the time evolution of the oc- 
cupation probabilities |cio| 2 and |coi| 2 , for different val- 
ues of g. In Fig. [21 we plot the results for very weak 
coupling, g/Ae = 0.03. The evolution takes the junction 
qubit 1 1) j j and transfers it to and from the resonator 
periodically. The exact, RWA, and dressed-state pertur- 
bative results arc all the same to within the thickness of 
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FIG. 4: Probabilities |cio| 2 and |cni | 2 for the strongly coupled 
case g/Ae = 0.30. Here both the RWA and dressed-state 
perturbative approximations fail. 



the lines shown in Fig. Thus, for this value of g, the 
RWA is extremely accurate. 

In Fig. 01 wc plot the probabilities for stronger cou- 
pling, g/Ae = 0.30. For this coupling strength, the 
RWA is observed to fail. For example, the RWA pre- 
dicts a perfect state transfer between the junction and the 
resonator, and does not exhibit the oscillations present 
in the exact solution. The dressed-state perturbative 
approximation docs correctly capture these oscillations. 
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FIG. 5: The fidelity of a state transfer from Josephson junc- 
tion to nanomechanical resonator as a function of interaction 
strength g. Note the difference in scale between the upper 
and lower curves. 

In Fig. 01 we show the same quantities for the case 
g/Ae = 0.5. At this coupling strength, both the RWA 
and the dressed-state perturbative approximation break 
down. 



IV. STATE TRANSFER FIDELITY 

In this final section, we briefly investigate to what ex- 
tent we may increase the junction-resonator coupling g, 
and still have an accurate state transfer from the Joseph- 
son junction to the resonator. As before, we start at time 
t — in the state |10). In order to define the fidelity of 
the state transfer operation, we first determine the time 
imin of the minimum of the probability |cio(i)| 2 Recall 
that |cio| 2 is the probability that the junction is in the 
m = 1 excited qubit state and the resonator is in the n — 
vacuum state. 

It will be convenient to define two fidelities: Fjj = 
1 — |cio(i m in)| 2 is the fidelity (or, more precisely, the fi- 
delity squared) for the junction, and F Tes = |coi(i m in)| ^ s 
the squared fidelity for the resonator S These quantities 



are different because of leakage to other states; however, 
in the RWA limit, they are both equal to unity. Fjj 
measure the success of de-exciting the qubit, and F rcs 
measures the success of exciting the resonator. In Fig. [3] 
we plot Fjj and F les as a function of g. Typically, the 
junction fidelity Fjj remains close to unity, with some 
oscillations, for all couplings. This behavior is a conse- 
quence of the fact that there is always a time where |cio| 2 
becomes small, as is evident in Figs. [3] and 01 However, 
because of leakage to other states, the resonator fidelity 
F rcs decreases significantly (again with oscillations due 
to the "switching" of i m j n with g) with increasing inter- 
action strength. The lower curve in Fig. [S] shows that 
F rC s > 90% is possible with g = 0.15 Ae, which allows a 
state transfer in under 5 ns. 

V. DISCUSSION 

We have developed a theoretical technique to ana- 
lytically calculate the leading-order perturbative correc- 
tions to the RWA or Jaynes-Cummings Hamiltonian for a 
quantum two-level system linearly coupled to a harmonic 
oscillator or single-mode boson field, a model central to 
many current quantum computing architectures. Such 
corrections are necessary to treat the fast information- 
processing regime where the interaction strength ap- 
proaches the qubit level spacing. The method was ap- 
plied to a current-biased Josephson junction coupled to 
a piezoelectric nanoelectromechanical disk resonator, and 
good agreement with exact numerical results was ob- 
tained. 
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If the higher lying junction states m = 2, 3, . . . are not ini- 
tially occupied, and the resonator also starts in the ground 
state, only the qubit states m = 0, 1 are coupled by "en- 
ergy conserving" terms. Therefore, the higher lying junc- 
tion states m = 2, 3, . . . can be treated with ordinary per- 
turbation theory in g. 

An alternative representation for the propagator is 

Gj/(t) = e - iw i t/h (^\Te- {l/h) f o dTV(T) \^,'), 

where V(t) = e iHjct/ft Ve- iH " ct/h , and where T is the time- 
ordering operator. However, directly computing G°y (t) by 
a series or cumulant expansion in V leads to a result that 
is valid only for short times. 

Because the desired final state is 1 01) , overall squared fi- 
delity of the state transfer operation is simply F TCE . 



